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Since Part II. of the Memoir appeared in November, 1898, the following papers by 
the author, bearing upon the Partition of Numbers, have been published : — 

" Partitions of Numbers whose Graphs possess Symmetry," ' Cambridge Phil. 

Trans./ vol. XVII., Part II. ; 
" Application of the Partition Analysis to the Study of the Properties of any 

System of Consecutive Integers," ' Cambridge Phil. Trans. / vol. XVIII. ; 
" The Diophantine Inequality Xx^: \iy" ' Cambridge Phil. Trans./ vol. XIX. ; 
"Combinatorial Analysis. The Foundations of a New Theory/' 'Phil. Trans. 

Roy. Soc. London/ A, vol. 194, 1900. 

In the present Part III. I consider problems of " Arithmetic of Position." In 
particular, I define a " general magic square " composed of integers and show that for a 
given order of square it is possible to construct a syzygetic theory. Such a theory is 
worked out in detail for the order 3 as an illustration. I further discuss the problem 
of the enumeration of the squares of given order associated with a given sum. I show 
that there is no difficulty in constructing a generating function for such squares even 
when the construction is specified in detail, and I obtain an analytical expression for 
the number when the sum, associated with rows, columns and diagonals, is unity 
or two. 

Art. 124. A " general magic square" I take to consist of n 2 integers arranged in a 
square in such wise that the rows, columns and diagonals contain partitions of the 
same number, zero and repetitions of the same integer being permissible among the 
integers. 

An ordinary magic square I define to be a general magic square in which the n 2 
integers are restricted to be the first n 2 integers of the natural succession. 

We may regard general magic squares as numerical magnitudes. To add two such 
magnitudes we add together the numbers in corresponding positions to form a 
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magnitude which is obviously also a general magic square. We can, therefore, form 
a linear function of magnitudes of the same order, n, the coefficients being positive 
integers, and such linear functions will denote a general magic square. 

The magnitudes, of the same given order, can be taken as the elements of a linear 
algebra, and since arithmetical addition can be made to depend upon algebraical 
multiplication, the properties of the magnitudes can be investigated by means of a 
non-linear algebra. 

Art. 125. The properties of a general magic square can be exhibited by means of 
homogeneous linear Diophantine equations, and it thence immediately follows that 
there must be a syzygetic theory of such formations. There exists a finite number of 
ground forms, corresponding to the ground solutions of the equations, and the method 
of investigation determines these and the syzygies which connect them. 

Generally speaking, there is a syzygetic theory associated with every system of 
linear homogeneous Diophantine equalities or inequalities, and it is because invariant 
theories depend upon such systems that they are connected with syzygetic theories. 

Art. 126. The method of investigation about to be given applies not only to magic 
squares of different kinds but to all arrangements of integers, which are defined by 
homogeneous linear Diophantine equalities or inequalities, whose properties persist 
after addition of corresponding numbers. 

For example, the partitions of all numbers into n 9 or fewer parts, are defined by the 
linear homogeneous Diophantine inequalities 

CCi — 0C 2 — «3 . . . — . Gfc w , 

and if another solution be 

Pi — H2 — Pz • • • — Hm 

we have 

and since the property persists after addition, a syzygetic theory results. 

This is one of the simplest cases that could be adduced and is at the same time the 
true basis of the Theory of Partitions. 

Many instances of configurations of integers in piano or in solido will occur to the 
mind as having been subjects of contemplation by mathematicians and others from 
the earliest times. These when defined by properties which persist after addition of 
corresponding parts fall under the present theory. 

Art. 127. There is no general magic square of the order 2 except the trivial case 

a ( it 

, but we may consider squares of order 2 in which the row and column 

properties, but not the diagonal properties, are in evidence. 
Let such a square be 

a x oc 2 
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which must clearly have the form 



a, a 



2 



oLo a. 



and we may associate with it the Diophantine equation 

«i + a 2 = a 5 

and regard a u a 2 and a 5 as the unknowns. 

The syzygetic theory is obtained by forming the sum 

for all solutions of the equation, and the result is 



a 



1 — 0/J\.i . 1 — CLjL-2 * 1 ^-5* 



where a is an auxiliary quantity and the meaning of the prefixed symbol O is that 

after expansion of the algebraic fraction in ascending powers of X 1? X 2 , X 5 we are to 
retain those terms only which are free from a. 
The expression clearly has the value 



1 



1 — -A.j A_5 . 1 — .A.2A.5 



The denominator factors denote the ground solutions 



«1 


a 3 


«5 


1 





1 





1 


1 



and the absence of numerator terms shows that there are no syzygies. 
Thus the fundamental squares are 



1 
1 



1 

1 



and this is otherwise evident. The case is trivial and is introduced only for the 
orderly presentation of the subject. 

Art. 128. Passing on to the general magic squares of order 3 we have the square 



a, a a< 



«4 % a 6 



CLrj OLq CCq 



40 



MAJOE P. A, MacMAHON: MEMOIR ON THE 



defined by the eight Diophantine equations 



ot 7 + a 8 + a 9 



= a 



10? 



a 



10) 



a 



10? 



= a 



10? 



a 2 + a 5+ a 8 

a 3 + a 6 + a 9 
ot 3 + % + a 7 



= a 



10? 



a 



10? 



a 



10? 



= a 



io» 



We require all values of the quantities a which satisfy these equations. 
To form the sum 

SX 1 ai X 2 a2 X 3 aa X 4 a4 X 5 a5 X 6 a6 X 7 a7 X 8 ^X 9 a9 X 10 a % 

for all solutions, introduce the auxiliary quantities 

a, 6, c, d, e, f, g, h 

in association with the successive Diophantine equations. The sum in question may 

be written 

1 



a 



r 



< 



( 1 - adgXJ ( 1 - aeX 2 ) ( 1 - afhX 3 ) ( 1 - bdXJ ( 1 - beghXJ 

X,„ \ 



(1 -6/Xe) (1 -c^X 7 ) (1 -ceX 8 ) (1 -cfgX 9 ) 1 - 



>• 



abcdefghj 



-j 



where after expansion we retain that portion only which is free from the auxiliaries. 
Remarking that 

1 1 



ft 



(1 -aPx) (1 -aP 2 ) (1 -aP 8 ) ( 1 - ^ 



PA - (1-P 1 P 4 )(1-P 8 P 4 ).(1-P8P4)' 



we eliminate the auxiliar a and obtain 



fl 



1 



1- 



X X X 



10 



-< 



bcefh, 



bcdfgh) \ bcdegj 



> 



( 1 - %AX 5 ) ( 1 - 6/Xe) ( 1 - cc^X 7 ) ( 1 - ceX 8 ) ( 1 - cfgX*) 



-j 



Put now bd = A, be = B, 6/' = O, cd — D, and we obtain 



ft 



r 



AXjX 



10 



BCD/* 



1- 



2^-10 

ODgh. 



X a X 



3-^MO 



BDy 



(1-AX 4 ) 



V 



< 



( 1 - Bc//iX 5 ) ( 1 - CX 6 ) ( 1 - D/*X 7 ) ( 1 - 



v. 



BDX 8 
A 



ODgX 

A 



an artifice which reduces the number of auxiliars to be eliminated by unity. 
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Remarking that 



(l-oPOd-aP^fl-ip. )(l->P. 



J. ~~" X jX 2-JL 3X 4 



1 PaP, 






(1 -P X P 4 ) (1 -P 2 P 8 ) (1 -P 2 P 4 ) ^ (1 -PxPs) (1 -P1P4) (1 -P 2 P 3 ) ' 



We eliminate A and find 



a 



/■ 



/ 






■< 



1 



:~jy- X-iXgXio 1(1 — JDJjXijLXg) (1 — 013^X4X9) ( 1 - 

~— — . X3X10 HI — IjffhjLfr j yi — OXg) (1 — JL//1X7) 



GDgh 



20x10 



> 



1 

J 



1 



jl n 



c& 










r / 1 \ f a \ i 1 \ 1 

[ - — -7^77 * OXiOXgOXio 1 ( X ^pry OX1OX9OX10 ! ( 1 Jl5 1^0X40x8/ (1 /^TY — T •^ k "2-^-10 I ! 



J 



V 



1 \ 

j^ — «— •"— — 0x30x10 ) \ *• ~*~ x~)Q flJSL.5 ! ( 1 v> oxg I ( 1 xJibJjL.'j} 



!► 



J 



Eliminating B from the first fraction and C from the second, we have 



O 



r 



i 



h 



^ 



< 



v 



v , rv7y \ I « "V "V 1 Y "V "V 

1 £ — Ijl/OAiAg) I 1 — UlvJS^iJ 1 x — - •— ■=—- — — - 0x20x10 j \ x irx 0x30V 50XIQ 

( y-sj-rT- \ / ^ 2"\7" "XT" XT" "17* \ f ~t -*- "XT* "XT'' "ct" "\7" 

x "**~ vy-A.g/ ( X """"* Q oV]-A-5xX9j\.j{) J I 1 "~" — • 0X30X40xg0X] 
\ g 



> 



4~ O 



gD 



Xv v y v 
I n * . '""'^ -mm Jfcw yi . Ji n . --IIn. V< i.4i. Jt» I I 9 



/• 



1 



^U V Y Y Y Y 1/1 f^Tl^^T Y \/l TiT) Y \ /l v xr 

0x10x40x^0x90x10 I \ 1 -""* \J.LJ(J 0x40x9 J ( x **~~ XJ/l/JjL.'j ) ! 1 —~ ">>j-r-^ — 7 0x20x10 



A 



A 



CD^/i 



> 



I l 1 "~" Vooxg J \ X "~~~ y oxioy^oxQoxjo I \ X oxgoxioxgox 

I ^ \ 9 



j 



+11 



h 



X 1 X 6 X 8 X 10 



r 



< 



l\ T3T iY Y \ /I TWo Y \/l a y Y Y Wl 

l X XJXJ 0x40X8 / \ X XJ(//t'Oxg I I X "-""" ^r~ "T" 0X2-^X5-^40 I 1 x 

\ DC/lb J \ 

j_ „^^ ^^ \ / 1 T~\ lb "ST \ i 1 "y "X7" "V "V \ 

=r — 0x30x10 J ( 1 XJfbj\.q\ ( 1 -—- — oxiOxgoxgoxiQ I 



^ Y Y Y 

cry- oxj^gox 



~\ 



10 



> 



EDg 



j 
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From the first fraction eliminate C, from the second C, and from the third B, 
obtaining 

1 



n 



r 



( 1 — DAX 7 ) ( 1 — jt-t X 2 X 6 X 10 j f 1 — =- X 3 X 5 X 10 



■\ 



> 



V 



(1 — g^XxXsXgXxo) f 1 — X 3 X 4 X 8 X 10 j f 1 — j X 2 X 4 X 9 X 10 



J 



+o 



gD 

h 



jLiJX.^JX.qjL.qJX. 



10 



r 



(l 



l 



DAX 7 )(l-^X 1 X 4 X 8 X 9 X 1 „j^l 



AoX«X 



2^6^M0 



/ 



< 



^ 



( 1 — Q XXX5A-9X10) i 1 X 3 X 4 X 8 X 

\ 9 



r 



1 — y X 2 X 4 X 9 X 10 
h 



/ ) 



y-XxXgXgXxo (1— g r X 1 X 3 X 4 X 5 X 8 X9Xio ) 



+0 



h 



r 



(1 -DKK,) ( 1 - ^X^X^xJ (l 



h 



Bgh 



XoXaX 



2-^6-*M0 



< 



\. 



jz X 3 X 5 X 10 ) (1 — <7 2 X 1 X 5 X 9 X 10 ) (1 — X 3 X 4 X 8 X 10 ) ( 1 



"> 



v. 



x "V "Y" "V" 



J 



Eliminating D from each of the three fractions, we obtain 



n 



r 



(1 — g^XiXgX^Xio) (l — X 2 X 6 X 7 X 10 j ( 1 — X 3 X 4 X 8 X 10 



< 



K. 



(1 — /^ 2 X 3 X 5 X 7 X 10 ) ( 1 — j X 2 X 4 X 9 X 10 



> 



J 



+ii 



1 



Xv "V "V "V "V V 
1 J\.2 -A 4 JSl. 6-A.g .A.9.A. 



10 



/" 



(1 — ^XxXsXgXio) ( 1 — -X 2 X 6 X 7 X 10 



-i W Y Y Y 

1 — — AoAiAsA 



^ 



9 



, 3 iL 4 iv 8 ii 10 



> 



V. 



1 

h 2 



y X 2 X 4 X9X 10 ) ( 1 — 7-5 X 1 X 2 X 4 Xg A.8 A9 A10 



^ 



-fft 



T XiXgXgXxo ( 1 — X 1 X 2 X 3 X 4 X 5 X 6 X 7 X 8 X 9 X 10 ) 
h 



r 



(1 — i 9 a X 1 X 6 X 9 Xi ) ( 1 — X 2 X 6 X 7 X 10 W 1 — X 3 X 4 X 8 X 10 



i 



V 



"I 



y 



(1 - A 2 X 3 X 5 X 7 X 10 ) (l - iXxXgXsXxo) (l - 1 XxX.X.XgXsXeXxo 2 ) 



J 
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Art. 129. Before proceeding to eliminate g and h, observe that if we now put 
g = h = 1 , we obtain the generating function for the solutions of the first six of the 
Diophantine equations corresponding to the squares which possess row and column 
but not diagonal properties. 

Putting g = h == 1, the generating function reduces to 



1 — .Ai_A-2-A-3-A-4-A_5 A.(5_A-7_A-8iY.9-A_ 



3 



10 



(1 — XiXsXgXxo) (1 — XxXeXgXxo) (1 — X 2 X 4 X9X 10 ) 
(1 — X2X6X7X10) (1 — X3X4X8.X40 ) (1 — X3X.5X.7XK)) 



indicating ground forms 



X1X5X9X10, 
XiXeXgXio, 

X 2 X4X 9 X 10 , 



X 2 X 6 X 7 X 10 , 

XY y Y 

X 3 X 5 X 7 X 10 



connected by the ground syzygy 



I-A.5-A.9-A.K) . -A^-A-g-A^-A-io • ^3^4-A.8-^VL0 — -A-i-A.g-A_g-A.xo . _AL.2-A_4.A-9.A_10 . -A.3_A_5_A-7_A.10j 

corresponding to the fundamental squares 



1 











1 











1 


1 














1 





1 









1 











1 


1 











1 





1 














1 



connected by the fundamental syzygy 



10 
10 
1 

10 
1 
10 



each side being equal to 



+ 



+ 



10 
1 
10 

10 
1 
10 



1 1 1 
111 
111 



+ 



+ 









1 


1 





.0 





1 











1 





1 





1 














1 


1 











1 











1 





1 





1 









This is the complete syzygetic theory of these particular squares of order 3. 

G 2 
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Art. 130. Resuming the discussion, we proceed to eliminate g and h and remark 
that the second fraction may be omitted as contributing no term free from h. 
Eliminating g from the first and g from the third, we have 



O 



1 "I" Jv-i J\.2&3&4&5&Q&y&$&9-&-iQ 



r 



/ 1 IT V 2V 2V IT 2V Y 3\ / 1 Y Y 2 V Y 2V 2V Y 3\ 

^1 — AiA.3 A 4 A 5 A 8 A 9 A 10 )\i — A.jA-2 A 6 A 6 A 7 A 9 A 10 j 



1 



(1 — A 2 X 3 X 5 X 7 X 1 o) f 1 — t X 2 X 4 X 9 X 10 



> 



v_ 



+n 



Y Y Y Y (^ Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 6 \ 

y- .A-j-A-gxVgxvjQ ( L -^-1 -A-2 -^-3 -^*-4 -A-5 -^t> -A-j -A-g -^-9 -^-VlO / 



/l Y Y 2 Y 2V Y 2V Y 3\ /I Y Y 2 Y Y 2Y 2V Y 3\ 

^ 1 — A x A 3 A 4 A 5 A 8 A 9 A 10 j [ 1 — A x A 2 A 5 A 6 A 7 A 9 A 10 j 



"^ 



< 



(1 — fb JL3A5X7.A40) ( 1 — y A1A.gA.3A. 



h 



I YYYYYYY 2 

•10 /l L T2 ^1^2-^4^-6^-8^-9^-10 



V 



Art. 131. If the diagonal property associated with g is alone to be satisfied in 
addition to the row and column properties we may put h = 1. Observe that the 
second of the three fractions cannot now be omitted. Simplifying we obtain 



1 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 

1 ^\-i -A_2 -^3 -^4 -A-5 -^6 -^7 -^8 -^-9 -^ 



10 



(1 — XiX 6 X 8 X 10 ) (1 — X 2 X 4 X 9 X 10 ) (1 — X3X 5 X 7 X 10 ) 

/I Y Y 2 Y Y 2 Y 2 Y Y 3\ / 1 Y Y 2 Y 2 Y Y 2 Y Y 3\ 

yl — Ai A.2 A5 Ag A 7 A. 9 Axo / \ 1 — A^ A3 A 4 A5A8 A9A10 J 



Establishing the five ground products 



AlAgAgAjo, 

a 2 A. 4 A. 9 A.i , 

A-3A 5 A 7 Ai , 

XY 2 Y Y 2 Y 2 Y Y 3 
\J\-2 -A.5^\.g Jx.'j J\.Q.j\-iQ > 

Xv 2 Y 2 Y Y 2 Y Y 



10 



connected by the ground syzygy 

(AiAeAgAioj (A 2 A 4 A 9 A 10/ | (A 3 A 5 A 7 A 10 j 

/Y Y 2 Y Y 2 Y 2 Y Y 3\ /v "V" 2 Y 2 Y Y 2 Y Y 3\ 

— ^AiA.2 A.5^\-g A 7 ^X-g^V-xQ J ^AjA3 A 4 -ZV5.A-3 -A. 9 ^V. 1 o J 

corresponding to the fundamental squares 



10 
1 
10 



1 

1 
1 



1 
1 
10 



12 
12 
2 1 



1 2 

2 10 
2 1 



THEORY OF THE PARTITIONS OF NUMBERS. 45 

connected by the fundamental syzygy 

100 10 001 
'2 001+2 100+2 010 
010 001 100 
12 10 2 
= 012 + 210, 
2 1 2 1 

involving the complete theory of the squares in which the property of one chosen 
diagonal is excluded. 

Art. 132. Resuming and finally eliminating h 9 we obtain 



1 + XiX.2-X.3X4-X-5X5X7X.gX9X; 



3 



10 



/l Y Y 2 Y Y 2 Y 2 Y "Y 3\ /1 Y Y 2V 2V Y 2V "V" 3\ /1 v 2y y 2 Y Y Y 2V 3\ 

[ 1 —^ViA.2 A 5 A 6 A7 A 9 A 10 j V 1 ""^1^3 A I A 5 A 8 A 9 A 10 ) \ L — A 2 A3A4 A 5 A 7 A 9 AlO ) 

X 2 Y Y Y 2 Y Y 2 Y 3 /l Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 6 \ 
1 u\.gJ\.K > J\_§ J\.>jJ\.g -^10 \ -*- — -^l -^*-2 -^-3 -^-4 -^5 -^-6 -^^7 -^-S -"^9 -^lO / 

(1 — X 1 X 2 2 X5X 6 2 X7 2 X 9 X 10 3 ) ( 1 — X 1 X 3 2 X4 2 X5X 8 2 X 9 X 1U 3 ) 

/I Y 2 Y Y Y 2 Y Y 2 Y 3\/i YYYYYYYYYY 3\ 

^1 — A x A3A5A6 -A.7A.8 A.10 ) ( 1 — ^iA.2-A.3A.4-A.5^6^7^g.A.9.A_i() J 



which may be written 



/I Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2 Y 2V ti\ 2 
[ 1 — A x A 2 A 3 A 4 A 5 A 6 A 7 A. 8 A 9 A 10 ) 



/I Y Y 2 Y Y 2 Y 2 Y Y 3\ /1 Y Y 2 Y 2 Y Y 2 Y Y 3\ / 1 y 2V Y Y 2 Y Y 2 Y 3\ ~> ' 

^1 -AiAg A 5 A 6 A. 7 A 9 A 10 ) [I — A!A 3 A 4 A 5 A 8 A 9 A 10 j ^1 — A x A 3 A 5 A 6 A 7 A 8 A 10 ) 

(1 — X 2 X3X4 X 5 X 7 X9 2 X 10 ) (1 — XXX2X3X4X5X6X7X8X9XJ0 ) 

indicating the ground products 

XY 2 Y Y 2 Y 2 Y Y 
lA 2 AgAg A7 A.9A] 

Xv 2V 2V Y 2 Y Y 3 
1A3 A 4 A 5 A 8 A 9 A 10 , 

X 2 Y Y Y 2 Y Y 2 Y 3 
l.A 3 A 5 A 6 A 7 A 8 A 10 , 



3 

-10 5 



X2V Y 2 Y Y Y 2 Y 3 
2 A3A4 A5A7A9 A 10 , 

XYYYYYYYYY 3 
x A 2 A 3 A 4 A 5 A 6 A 7 A 8 A 9 A 10 

connected by the fundamental sy zy gies 

/Y Y 2 Y Y 2 Y 2V Y 3\ /y Y 2 Y 2 Y Y 2 Y Y 3\ 
\A X A 2 A 5 A 6 A 7 A 9 A 10 ) ( A X A 3 A 4 A 5 A 8 A 9 A 10 / 

— /Y 2 Y Y Y 2 Y Y 2 Y 3\ /y 2V Y 2 Y Y Y 2 Y 3\ 

— {A x A 3 A 5 A 6 A 7 A 8 A 10 ) (A 2 A 3 A 4 A5A7A9 A 10 / 

(y YYYYYYYVY3\2 

— ^A 1 A 2 A3A4A 5 A 6 A7A 8 A9^V 10 J 

corresponding to the fundamental general magic squares 
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2 1 2 1 

12 2 10 

12 10 2 



1 


2 





1 





2 









1 


2 


2 


1 









2 





1 





2 


1 
1 
1 
1 


1 
1 
1 


1 
1 
1 



connected by the fundamental syzygies 

120102 111 201 21 
012 + 2 10=2111= 012 + 210. 
201 021 111 120 102 

Art. 133. If the sum of each row, column, and diagonal be Sn, the number of 
general magic squares of order 3 that can be constructed is, from the generating 
function, the coefficient of a? 3M in the expansion of 

and this is found to be 

n 2 +(n+lf. 

Art. 134. The ordinary magic squares, the component integers being 0, 1, 2, 3, 4, 
5, 6, 7, 8, are eight in number and are easily found to be 

120 201 723 

012+3 012 = 048, 
2 1 12 5 6 1 

and seven others obtained from 

120 201 120 021 120 02 1 
3012 + 012, 012+3 210, 3 012 + 210, 
201 120 201 102 201 102 

102 201 102 201 102 021 
2 10+3012, 3210 + 012, -210+3210, 
021 120 021 120 021 102 

2 1 2 1 
3210 + 210. 
10 2 10 2 

Art. 135. There is no theoretical difficulty in proceeding to investigate the squares 
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of higher orders, but even in the case of order 4 there is practical difficulty in 
handling the H generating function. There are 20 fundamental squares, viz. : — 



10 
10 
1 
10 

10 

10 

1 

1 

110 
10 10 
10 1 
11 

2 

10 10 

11 

10 1 

10 1 

110 

10 1 

2 



10 
1 
10 
10 

10 
10 
1 
10 

10 10 
11 
110 
10 1 

2 

10 1 

110 

10 1 



10 1 

11 

10 10 

2 



10 
10 
10 
1 

1 
10 
10 
10 

10 1 
110 
11 
10 10 

10 10 
2 
110 
110 

11 

110 

2 

10 1 



10 
1 
10 
10 

1 
10 
10 
10 

11 
10 1 
10 10 
110 

110 
110 
2 
10 10 

10 1 

2 

110 

11 



§10. 

Art. 136. The direct enumeration of general magic squares of given order and sum 
of row. 

Let h w denote the sum of all the homogeneous products w together of the 
magnitudes 

If h w be raised to the power n and developed, the coefficient of 

is the number of squares that can be formed of order n 9 so that the sum of each row 
and column is w 9 but in which there is no diagonal property in evidence.*" 

* "Combinatorial Analysis—The Foundations of a New Theory," * Phil. Trans./ A, vol. 194, 1900, 
p. 369 et seq. 
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In fact, if 

(x— ol x ) (x — a 2 ) t ..(x—a n ) = af t —p l a?~ 1 + ... +( — ) n p n 
and 

w I D w = (d Pl +p x d p2 + . . . +p n ^ d pn f\ 

an operator of order w obtained by raising the linear operator to the power w 
symbolically as in Taylor's theorem, then the number in question is concisely 
expressed by the formula 

a particular case of a general formula given by the author (loc. cit.). 
Art. 137. To introduce the diagonal properties, proceed as follows: — 



Let hj s) denote what h w becomes when Xa s , jx,a n __ s+1 are written for a s) a n _ s+1 

respectively, and form the product h w a) hj 2) . . . hj n \ 

I say that the coefficient of 

(\fji(X 1 a 2 ...a n ) w 

in the development of this product is the number of general magic squares of order n 
corresponding to the sum w. 

To see how this is take n = 4, w = 1, and form a product 

x (aj + Xa 2 + /xa 3 + a 4 ) 
x (a x + /xa 2 + Xa 3 + a 4 ) 
x (/xa x + a 2 + a 8 + Xa 4 ), 

and observe that, in picking out the terms 

one factor must be taken every time from each row, column and diagonal ot the 
matrix. 

Similarly, if n = 2, we form the product 

{ X 2 «! 2 + X/xa^ + /x 2 a 4 2 + (Xa x + jua 4 ) (a 2 + a 8 ) + a 2 2 + a 2 a 8 + a 3 2 } 
x {X 2 a 2 2 + XjLta 2 a 3 + ju 2 a3 2 + (Xa 2 +/xa 3 ) (a 1 + a 4 ) + a i 2 + ^i a 4 + a 4 2 } 
x {X 2 a 3 2 + X/xa 2 a 3 + /x 2 a 2 2 + (Xa 3 +/xa 2 ) (a x + a 4 ) + a x 2 -f a^H- « 4 2 } 
X {X 2 a 4 2 + Xjiia 1 a 4 + /x 2 ai 2 +(Xa 4 + />ta 1 ) (a 2 + a 3 ) + a/ + a 2 a 3 + a 3 2 }. 

In forming the term involving 

X 2 /x 2 a 1 2 a 2 2 a 3 2 a 4 2 

regard the successive products as corresponding to the successive rows of the square, 
the suffix of the a as denoting the column, and X, fx as corresponding to the 
diagonals. 

Thus picking out the factors 

X 2 a x 2 5 /xa 3 a 4 , /^a 2 a 4 , a 2 a d , 
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we obtain the corresponding square 

2 

11 

10 1* 

110 

These examples are sufficient to establish the validity of the theorem. 
Art. 138. If we wish to make any restriction in regard to the numbers that appear 
in the s th row, we have merely to strike out certain terms from the function 

E.g., if no number is to exceed t, we have merely to strike out all terms involving 
exponents which exceed t. 

If the rows are to be drawn from certain specified partitions of w 9 we have merely 
to strike out from the functions 

h CO l (2) J, (n) 

all terms whose exponents do not involve these partitions. 

We have thus unlimited scope for particularising and specially defining the squares 
to be enumerated. 

Let us now consider the enumeration of the fundamental squares of order n 9 such 
that the sum of each row, column and diagonal is unity. Observe that if the 
diagonal properties are not essential the number is obviously n ! 

Art. 139. It is convenient to consider a more general problem and then to deduce 
what we require at the moment as a particular case. I propose to determine the 
number of squares of given order which have one unit in each row and in each column, 
and specified numbers of units in the two diagonals. 

Consider an even order 2n, and form the product 

x (a x + Xa 2 + Og + . . . + a 2w _ 2 + /xa 2w _ 1 + ot 2n ) 
x (fl^ + a 2 + Xa 3 + . . . + /xa 2w _ 2 + a 2n -i + a 2n) 



x (a x + a 2 + fia s + . . . 4- Xa 2w _ 2 + a 2n _ x + a 2n) 

x (ax + fia 2 + a 3 + . . . + ot 2n _ 2 + Xa 2n _! + a 2n ) 

x (/xa-L + a 2 + a 8 + . . . + ot 2n _ 2 + a 2w _ L + X.a 2w ). 
We require the complete coefficient of 

when the multiplication has been performed. 
Writing ^ 

VOL. CCV. — A. H 



50 MA JOE P. A. MacMAHON: MEMOIE ON THE 

the product is, taking the t th and 2n+l— t th factors together, 

{s + (X-l)a 1 + (ft-l)a 2w } {5 + (/Lt-l)a 1 + (X-l)a 2w } 
x{5 + (X~l)a 2 +( i u,--l)a 2?l _ 1 } {5 + (^-l)a 2 + (X-l)a 2w _ 1 } 
x{s+(\-l)a 3 +(/x-l)a 2M „ 2 } {5 + (/x-l)a 8 + (X-l)a an . 2 } 

x{5 + (X— l)a„+(/x— l)o„ +1 } {5 + (ft— l)a B + (X— l)a B+1 }. 

Observing that we only require terms which involve the quantities a with unit 
exponents, the product of the first two complementary factors is effectively 

^ 2 + (X + /x-2)(a 1 + a 2n ).9+{(X^l) 2 +(/x-~l) 2 }a 1 a 2n , 
and the complete product 



t—n 

n 

4=1 



n[^+(X + /x-2)(^ + a 2 , +1 _,)5+{(X-l) 2 + (^-~l) 2 }^a 2w+1 „J 



has, on development, the form 

a** + A r 9 2w_1 + A 2 s 2w ~ 2 + . . . + A 2 „, 

where A w is a linear function of products of the quantities a, each term of which 
contains m different factors a, each with the exponent unity. 
Since, moreover, s m gives rise to the term 

ml Sa 1 a 2 ...a m , 

it follows that the coefficient of 2a 1 a a ...a 2n in the product is obtained by putting each 
quantity a equal to unity and s m = m !. 

Hence, if s w = m! symbolically, the symbolic expression of the coefficient is 

{s 2 +2(\+fx-2)s + (\-l) 2 + (iM-l) 2 y\ 

or 

{6« 2 -4.9 + 2 + 2(X + /x)(5-l) + X 2 + /x 2 } w , 
or writing 

s 2 — 4.9 + 2 = cr 2 , $—1 = cti 

{cr 2 + 2(X + /x)cr 1 + X 2 4-/>t 2 } w . 

This is the complete solution of the problem for an even order 2n. 
For an uneven order 2w+l, it is now evident that the symbolical expression of the 
coefficient of 

IS 

{cr 2 + 2 (X + /x) cri + X 2 + ^ 2 } w (a-i + X/x), 

the complete solution in respect of the uneven order 2n+l. 
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Art. 140. To find the number of ground " general magic squares" corresponding to 
the sum unity, we have merely to pick out the coefficient of X/x ; we thus find 

even order 2n number is 8 ( „ J cr 2 w ~ 2 cr 1 2 ) 

uneven order 2n+l number is 8 ( i? ) cr 2 H "" 2 cr 1 8 +cr 2 tt , 

wherein it must be remembered that the a products are to be expanded in powers 
of 6* and then s m put equal torn! 

In the general results the coefficient of 

xy i 

gives the number of squares in which the row and column sums are unity and the 
dexter and sinister diagonals' sums are I, m respectively. 
I give the following table of values of simple <x products : — 







&l 






01 





r 2 








g 


cricr 2 








Vl 


cricr 2 


&2 






O-l 5 


or*ar 2 


(r x cr 2 




o-i 6 


4 

o"! <r 2 


2 2 


*■ 


7 
<7i 


(Tx<T 2 


~ 3^. 2 
(Tx <T 2 


CTjCTa 3 














1 











2 











9 


4 








It—Mir <— lir 


* 


16 




265 


168 


116 


80 


1854 


1280 


920 


672 



The numbers erf = (s—l)*' denote the number of permutations of ^ letters in which 
each letter is displaced and constitute a well-known series. 

The remaining numbers are readily calculated from these by the formula 

d (T 2 — (T\ <T 2 — 6<Ti <J 2 — (T1CF2 . 

Art. 141. Another solution of the same problem yielding a more detailed result is 
now given. 

For the even order 2n I directly determine the coefficient of 

A JJL 0tj0t 2 • • * GL 2n 

in the product above set forth. 

H 2 
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We have to pick out I X's and m /x's and to find the associated factors, 2n—l—m in 
number, which are linear functions of the quantities a. 

In any such selection of I X's and m /x's there will be i pairs of X's symmetrical 
about the sinister diagonals and j pairs of /x's symmetrical about the dexter diagonals, 
and the associated factors will depend upon the numerical values of i and j. 

Consider then in the first place the number of ways of selecting I X's in such wise 
that i pairs are symmetrical about the sinister diagonals. 

This number is readily found to be 

i)\l—2i) 

With these I X's we cannot associate any /x which is either in the same row or in 
the same column as one of the selected X's. 

Each of the i symmetrical pairs of X's in this -way accounts for 2 /x's, and each of 
the 1—2% remaining X's accounts for 2 /x's. 

Thus we must select m /x's out of 2n-2i—2 (l—2i) /x's, i.e., m /x's out of 
2n—2l + 2i /x's. 

We may select these so as to involve i pairs symmetrical about the dexter 

diagonals in 

J A m-2j / J 

This number is obtained by writing in the first formula n—l~\-i,j and ra for 

n 9 i and I respectively, 

and observe that we may do this because the selection of a symmetrical pair of X's or 
of one of the remaining X's results in the rejection of a pair of /x's which is 
symmetrical about the dexter diagonals. 

Consequently the 2n~ 2l + 2i possible places for the m /x's are also symmetrically 
arranged about the dexter diagonal Hence the formula is valid. 

We have established at this point that we may pick out I X's involving i 
symmetrical pairs and m /x's involving j symmetrical pairs in 

(?) (riK* • ( B -j + *) (r^-ifr" ^ 

We must now determine the nature of the 2n—l—m associated factors, linear 
functions, of the quantities a. 

In the matrix of the product delete the rows and columns which contain selected 
X's and /x's. We thus delete l + m rows and l + m columns. 

Consider the 2n—l—m remaining rows. There remain in these rows at most 

2n~-l—m 
elements a, because l + m columns have been deleted, but some of these elements 
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must be rejected if they involve X or //, as coefficients, because by hypothesis we are 
only concerned with I X's and m /x's, and these have already been accounted for. 

Observe now that the columns which contain a symmetrical selected pair of X's 
only contain jit's which are in the same rows as these X's, and therefore the deletion of 
these columns cannot delete /x's appertaining to any rows except those occupied by 
the selected pair of X's. Observe further that the column which contains an 
unsymmetrical X, say in the p th row, contains a /x in the 2n—p + l th row, and that 
therefore the disappearance of a /x in the 2n— p+l th row follows from the deletion of 
a column containing an unsymmetrical X in the p th row. 

Hence of the 27i — l—m rows in question 

l + m—2i—2j rows contain 2n—l—m—l 9 a elements, 
and thence 

2n — 2l— 2m+2i + 2/ rows contain 2n~l—m—2, a elements. 

Accordingly if s is the sum of all the a elements except those which appear as 
coefficients of the selected X's and /x's the co-factor of 

n\/n-~i\ 2l - 2i /n-l + i\/n-l + i-j\ 2m -. 2j 
i)\l-2i) m \ j )\ m-2j ) 

contains I + m — 2% — 2j factors of type 
and 2n—2l—2m+2i+2j factors of type 

or of n— I— m + i+j squared factors of type 

(s-a v -a w )\ 

since these factors occur in equal pairs. 

Hence the co-factor is 

n (s - a u ) n (s - a, - a w ) 2 , 

wherein the quantities a M , l + m— 2i— 2j in number, which appear in the first product, 
and the quantities a V9 a, m 2n~2l—2m + 2i + 2j in number, which appear in the second 
product, are all different. 

Also (s— a v — ol w ) 2 is effectively equal to 

s 2 ~- 2 (a v + a w ) s + 2a v a. w 

since squares of the as may be rejected. 

Hence, by the reasoning employed in the first solution we may put the quantities 
a equal to unity, regard s p as equal to p ! symbolically, and say that the coefficient of 

in the product 

U(s-a^n{s 2 --2(a v + a w )s + 2a v a w } 



04J: 
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has the symbolical expression 

or, putting 

«s— 1 = <Ti, s 2 -~ &s + 2 = <r 2 , 
we obtain 



n\ /n-A 9 ^ 2i f?i~~l+i\ ( n -l + i-f\ m 



l-2i 



02 



J? 



m— -2? 



Oi-2j l*f m— 2i~2j n~l~~ni-¥i-¥j 



for the number of squares such that — 

(1) Sum associated with rows and columns is unity ; 

(2) There are I units involving i symmetrical pairs in the dexter diagonal ; 

(3) There are m units involving/ symmetrical pairs in the sinister diagonal 

Giving i and j all possible values we find that the complete coefficient of 



^1^2 * * * ^2; 



n 



in the product, which we have already ascertained to have the expression 

{o- 2 +2(\+/x)cr 1 + X 3 + / i 2 } B , 
may be also expressed in the form 



I m 



f) (f) 2 l . ( n ~ l ) ( n Z V )V*<ri* m <r2- l ~ m 



m 



+ ( ? ) (^^1^) 2 l "~ 2 J ^~i + "^ ( n ~~t+ 1 J g^ l+w-2_ »-*-»»+l 



l/W-2 







\ 



m 



<3*2 



™T* 



n\ n 







^19^ / ^ \( n ~~ M Om-2^ l+m~2^ n-Z-m+1 



+j 



(n\ fn—2\^i^^ In— 1 + 2^ ln—l + 2^ m ^i +m ^ i _ n ^ l ^ m+2 



\ 



2/W-4 







m 



<*2 






rA in—1 \ g^-s ln—l+ 1\ / n— -6 * 2 TO ~* 3 rr *+*»~<l_ n-*-?»+2 



1 / W-2 



m—2 



2 



4j (Tj 



+l 6;r/; 2 -v 2 



m—4t 



CTc 



+ (J) (r 6 3 ) 2l " * (* "o + 3X ) ^ + 3 ) 2^*— err '- +8 



2/U-4 



m—2 






0* f 



+ 



7&\ /% — 1 \ 2^—3 /^— ^+1\ M — £— 1 \2»~4 Z+m-6 w-Z-m+3 



l/W-2 






m— 4 



+ (!^ f ^ 2 1 , /^r^ /^~**~ 3 ^ 2 w "V 1 ' +w "Vr '~ w+8 



\0 \l 



3 



m— 6 



! » » • 



A. fJb » 
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Further simplification of this series cannot be effected because each term of the 
sum must be considered on its merits and does or does not add to the numerical 
result as may appear. 



Art. 142. Writing the result for even order 2n 



{<x 2 + 2 (\+ fi) o-! + X 2 + jU, 2 }" 

it appears that the result for uneven order 2n + 1 may be written 

{ a, + 2 (X + fi) o-! + X 2 + fx 2 } n (o-j + V) 

= t${<r$(n, l,m) + F(n, Z-l, ro-l)}xy\ 

For the squares of simple orders we have the results — 



II 



Order 2. 



1= 















1 
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Z= 
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Order 3. 



1 






2 









2 
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1= 
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s 



Order 4. 



1 
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4 
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Order 5. 






16 


16 
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4 


8 
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16 
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4 





2 
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Order 6. 



1= 



6 






80 


96 


60 


16 


12 

3 



0. 







1 







1 


96 
60 
16 


96 


48 

24 




24 
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48 
24 
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12 
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Art. 143. I now proceed to consider the enumeration of the squares of even 
order 2n, such that every row and column contains two units, and the dexter and 
sinister diagonals I and m units respectively. 

I form the product 



(1) (2) a (2n) 
^2 ^2 • • • ^2 ; 



where a 2 {s) is the sum two together of the quantities 

and I seek the coefficient, a function of X and /x, of 

y0t>i&> 2 * ..&> 2n ) 
in the product. 

The coefficient of X z ju w in the sought function of X and fi is the required number. 

Let p u p 2 be the sum and the sum two together of the quantities 



then 



a 



a) 



a 



(2n) 



whence 



a 2 (1) a 2 (2n) 



^l> ^2? ,,,a 2»5 

= ^2+{(^-l)a 1 + (\-l)a 2 4(Pi~ a i^ a 2n) + (X/x--l)a 1 a 2H , 

: p 2 2 +(k + fjL-2) (ot 1 + a 2n )p 2 p 1 -(\ + [i — 2) (Gt l + a 2n fp 2 
+ 2 (X/x, — 1 ) a x a 2n p 2 + (X/x — 1 f a x 2 a 2n 2 
+ (X-l)( / x^l)K + a2 /) + {(X-l) 2 + (/.-.l) 2 }a ia2riJ 9 1 2 
- 2 (X 2 + Xft + fi 2 - 3X - 3/x, + 3) { cL Y a 2n {a x + a 2?i ) p x - a^a^ 2 } 
+ (X//,-l) (X-f/ut— 2) a 1 a 2ll (a 1 + a 2ll )jp 1 -2 (Xft-l) (X + /A-2) a x 2 a 2? ; 
+ terms involving powers of a l5 a 2n above the second. 



.2 
2 
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The product of a 2 {1 \ a 2 {2n) is thus, after re-arrangement, effectively equivalent to 

p 2 2 + (X + /x - 2) ( ai + a 2w ) p 2 ^^ 

+ 2 (X-l) (ft — 1) aia Jte jp a + (X-l) (/a— 1) (<h 2 + a 2 ?)p{ 

+ {(\-l) 2 + {[i-l) 2 }ct 1 a 2)l p 1 2 

+ {(X + ^)(X/a-X--/a)-(X + ^ — 1)(X + /a— ^laiag^ai + aa^)^ 

+ {(X/x-X-/*) 2 + (X + /a- 1) (X +ft-3)} a^ag/. 

Regarded apart from jp 3J .Pi this expression is a function of a x , a 2/l ; the product 

a 2 ^a 2 (2 ^ l) 
is a function of a 2 , a 2?i -i, and generally the product 

is a function of a„ ot 2? , +1 _ 5 , and all of these products are of similar form in regard to 
Remembering that we desire the coefficients of 

in the product 

a (1) r/ (2) a & l) 

U/ 2 U/ 2 • • • U/ 2 y 

we must distinguish between p 2 where it occurs as a multiplier of a 2j ru 2 2 and where 
it occurs as a multiplier of a y a 2m and make a similar distinction in respect of p 2 . 
Put then 

(<*i+*2?)Pi = (ai 2 +«2» 2 ) nA 

Putting further the quantities a equal to unity and regarding a product 
as a symbol for the coefficient of symmetric function 

/oa+rf-i & + 2c\ 

in the development of symmetric function 

(l 2 ) a+c (l) 6+arf , 
I say that 

[p 2 2 + 2(X + /x~2)p 2 p 1 ~2(X + ^^2)|9 2 +2(X~-l)( / x-l)(7r 2 + 7r 1 2 ) 

+ {(X-l)V0*-l)'} i > 1 ' 
+ 2{(X + /x)(X/x-X- / x)-(X + /x-l)(X + /x,-4)}jp 1 

+ (X / £-X~/x) 2 + (X + /x-l) (X + /x-3)] ;i 

is the symbolic expression of the required coefficient of 

yCt l (X 2 . , ,w 2n y. 
VOL. CCV. A. I 
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This may be written 



where 



{o- 4 + 2 (\ + il) o- 3 +(X 2 +/* 2 ) cr 2 +2X/xc7 / 2 + 2X i a (X+/x) c^ + X 2 /* 2 }", 



0" 4 

CT 3 
<T 2 



p 2 2 ~4p 2 j9 1 + 4(p 2 +7r 2 ) + 2( : p 1 2 +7r 1 2 )~8^ 1 + 3 5 

P2P1- (i>2+^ 2 ) ~ (i>l 2 + 7T 1 2 )+ 5^-2, 

Pi 2 -4^ + 2, 

7T2 + 7T! 2 — 4^ + 2, 
_Pl— 1. 



For the uneven order 2n + l it is easy to show that the coefficient is symbolically 
{or 4 +2 (X + /x) o- 3 + (X 2 + ^ 2 ) cr 2 + 2X/tcr / 2 + 2X/x(X + ^)cr 1 + XV 2 } n x(^ 2 -cr 1 + a- 1 A/x). 



It is easy to calculate the values of 



p 2 a Pi^2^i d 



for small values of a, &, c ? d. 

Some results are, omitting the obvious result pf = & !, 
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enabling the. verification of the results 



0-4 = 0-3 = 0-2 = <t' 2 = cr x = 0. 



cr 2 2 = 4, 0-3^ = 0, a- 



/ 2 



% 



Hence for the even order 2 the whole coefficient is \ 2 /x 2 5 corresponding to the only 
possible square 



1 


1 


1 


1 



and I find for the uneven* order 8 

\ 2 +/?+2\ 2 //, 2 (\ + /u). 

Art. 144. To find in general the number of squares which have two units in each 
diagonal we find the coefficient of X 2 /x 2 and obtain for even order 2n. 



I) <r " 1+ [2) ^" 2 (2<r/ + 4o-?+ 16CT8O-0 
+ (£j a-r 3 (24cr3 2 (T 2 + 48cr 8 V a ) + (fy a-/"" 4 96o- 3 4 ; 



putting n = 2 we find for the order 4 



2cr 4 +2cr 2 + 4cr 2 + 16cr 3 (7i == 16. 

and the verification of this number is easy. 

For the uneven order 2n+l we obtain the number 



^ »-i »~* ~ -1- ( n \ <?*-*%<? fa 



1 o-± .^o" 2 cr 1 +( 2 



n\ n-L^/n 



+ (P2-<ri)l['{)<rr l + 2 K2 



o-/~ 2 (o- 2 2 +2o- / 2 2 +8a- 3 cr 1 ) 
(ri\ 



+ ( Jj a/- 3 24 (<rfa+2<rfa 2 ) + Q a-/"" 4 . 96a-, 



The general value of 



AW 



may be obtained by means of the calculus of finite differences. 

There is no theoretical difficulty in finding symbolical expressions for the enumera- 
tion of general magic squares associated with higher numbers, but the method does 
not lead to the determination of general magic squares. These must be regarded as 
arising from the generating function method of § 9. 
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